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Abstract 

Background: Hybridization and polyploidy are central processes in evolution and speciation. These mechanisms 
often lead to complex patterns of genetic variation and the creation of novel genotypes, which may establish if 
they become isolated from gene flow. However, in the absence of reproductive isolation, species boundaries might 
easily be disrupted. Here, we used a combination of AFLPs, chloroplast DNA markers and flow cytometry to 
investigate the evolutionary outcomes of hybridization between two endemic Ecuadorian species of Epidendrum 
(£ madsenii and £ rhopalostele) in three hybrid zones. Postmating isolation was also quantified to determine the 
role of this barrier in restraining gene flow between hybrids and the parental species. In addition, future ecological 
niche models were constructed to predict the outcomes of hybridization between these species. 

Results: Our results confirmed the presence of hybrids in all hybrid zones, but revealed that a third parental species 
(£ falcisepalum) has contributed to one of the hybrid zones studied. Backcross genotypes were frequent in all 
hybrid zones, which was in accordance with the absence of strong reproductive barriers. The process of 
hybridization was highly asymmetric and followed in some cases by polyploidy. The projection of future niche 
models predicted a severe reduction in the area suitable for the occurrence of these species, although favorable 
conditions will still occur for the existence of the current hybrid zones. 

Conclusions: The recurrent process of hybridization has compromised the genetic integrity of the parental species. 
Most individuals of the parental species can no longer be considered as pure-bred individuals because most were 
classified as backcrossed hybrids. Novel genetic lineages occur in all hybrid zones implying that hybrids are fertile 
and can compete with the parental species. These results, together with the prediction of suitable conditions for 
the future occurrence of these hybrid zones, highlight the importance of conserving these geographic areas as 
sources of novel taxonomic entities. 
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Background 

Natural hybridization has been reported to occur in a 
variety of organisms [1], and is often associated with 
areas where previously isolated lineages come into con- 
tact and mate, producing offspring of mixed ancestry 
[2], If hybrids are formed easily, a hybrid zone may de- 
velop, which may persist over time [3] or lead to the 
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fusion of parental species into a single interbreeding spe- 
cies [4,5]. 

Introgression with one or both parental species can 
ease the progress of genetic swamping [6,7], which 
coupled with hybrid heterosis [8,9] can further enhance 
the process of homogenization or displacement of paren- 
tal species. In some cases, introgressed genotypes colonize 
new habitats or become genetically stabilized, providing a 
pathway to the evolution of new lineages [10]. Conversely, 
strong isolation or strong selection against hybrids can fa- 
cilitate assortative mating and reinforce reproductive bar- 
riers between parental hybridizing species [11,12]. 
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The frequency and success of interspecific matings 
therefore influence the dynamics and evolutionary out- 
comes of hybridization. In some contact zones, intermedi- 
ate hybrid genotypes predominate ('unimodal hybrid 
zones'; [13,14]), whereas in others, hybrids are rare and 
populations consist largely of individuals genetically simi- 
lar to one or other parental genotype ('bimodal hybrid 
zones'; [15]). Nonetheless, studies in several different or- 
ganisms demonstrate that genotype frequencies may vary 
considerably within the same hybrid zone and create a 
continuum from unimodal to bimodal genotype frequen- 
cies [16,17]. Finally, species might geographically overlap 
but never hybridize, mostly because of strong prezygotic 
isolation [18]. 

Understanding the strength of reproductive barriers 
might provide important clues concerning the evolution 
of hybrid zones [19,20]. Prezygotic barriers are usually 
considered to be most important in the process of speci- 
ation [21-24]. Asymmetries in the flowering period of 
co-occurring species or specificity of pollinators are 
among the most important prezygotic barriers [25] and, 
indeed, might shape the structure of hybrid zones. For 
instance, bimodal hybrid zones are usually characterized 
by having strong prezygotic barriers, in contrast to their 
weakness in unimodal hybrid zones [17]. However, post- 
zygotic barriers might also influence the outcome of 
hybridization, especially in cases where prezygotic bar- 
riers are incomplete and allow frequent homogenizing 
gene flow [26,27]. 

Polyploid formation, which was formerly expected to 
promote immediate reproductive isolation between the 
new incipient species and its progenitors [22], is now- 
adays considered to be more complex than originally 
thought [28]. For instance, autotetraploids of Chamerion 
angustifolium are not instantly isolated from their dip- 
loid progenitors, although isolation might arise in time 
[29,30]. Furthermore, asymmetric gene flow from the 
diploid parent Capsella rubella to its allotetraploid pro- 
geny C. bursa-pastoris has been demonstrated and con- 
tributes significantly to high genetic variation in the 
novel polyploid [31]. Because a single origin of a poly- 
ploid represents an extreme bottleneck in this complex 
pattern of speciation, multiple scenarios might also be 
invoked to understand polymorphism across ploidy 
levels such as for instantaneous, recurrent polyploidiza- 
tion or reproductive isolation of polyploids [31]. 

Orchid species (Orchidaceae) are particularly prone to 
hybridization most likely because of the high number of 
sympatric species, in combination with a general lack of 
complete reproductive barriers [10]. Many studies, espe- 
cially those devoted to Mediterranean orchids, support 
this notion, and have provided many examples of how 
hybridization between orchids is influenced by several pre- 
zygotic [32,33] and postzygotic barriers [34]. In addition, 



vegetative spread and polyploidy may contribute to 
stabilization of many of these hybrids [35]. Outside the 
Mediterranean region, one of the best-known examples 
where hybridization is thought to have a strong influence 
in diversification is the genus Epidendrum. It is the largest 
orchid genus in the Neotropics and comprises about 1500 
species distributed from the United States to Argentina 
[36]. Plants of intermediate morphology occur in many 
populations and the existence of many hybrid zones have 
been reported based on morphological data and geograph- 
ical overlap between different species of Epidendrum 
[36-38]. Understanding the process of hybridization in 
Epidendrum based on morphological characters is a 
challenging task owing to the high variability of spe- 
cies, which often hinders reliable identification of spe- 
cies in the field [38]. However, despite the general 
belief that hybridization is frequent, only one molecu- 
lar study has reported the consequences of hybridization 
in Epidendrum [39]. 

In this study, we sought to understand the dynamics of 
gene flow between Epidendrum madsenii Hagsater and 
Dodson and E. rhopalostele Hagsater and Dodson, two 
epiphytic deceptive orchids that co-exist in three hybrid 
zones in southern Ecuador [36]. Epidendrum madsenii 
and E. rhopalostele belong to two different taxonomic 
groups and are well differentiated by their morphological 
traits [36]. The occurrence of hybridization between these 
species has been postulated on account of the existence of 
individuals with intermediate morphological traits in some 
areas of sympatry (Riofrio et al, unpub. data). Most spe- 
cies of Epidendrum have low pollinator specificity, which 
together with a high reproductive compatibility even 
between unrelated species suggests that hybridization 
may be common [39]. The main question addressed in 
this study was whether morphologically intermediate 
individuals were indeed hybrids between E. madsenii 
and E. rhopalostele, or simply fell within the wide range 
of morphological variation shown by these species. 
Given that other congeneric species are also present in 
the same area and these species are pollinated by gener- 
alist insects, the possibility that different parental spe- 
cies might give rise to morphologically similar hybrids 
must also be considered [40]. 

To answer these questions, we investigated the genetic 
structure of the three known hybrid zones and the nature 
of the morphologically intermediate specimens using an 
approach that combined AFLPs, cpDNA markers, genome 
size and controlled pollinations in conjunction with 
climate-based predictive models. The results allowed us to 
(1) characterize the genomic composition of hybrid zones 
and their hybrids (e.g., ¥ 1 hybrids and backcrosses); (2) de- 
termine the presence of asymmetrical hybridization pat- 
terns and the dynamics of genome size in different hybrid 
zones; (3) investigate the strength of reproductive barriers 
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in different hybrid zones; and (4) based on the molecular 
data and niche models, predict the long-term fate of the 
hybrid zones between these deceptive orchid species. 

Results 

Plastid DNA diversity and haplotype network 

Statistical parsimony analysis yielded a single network 
containing two major groups with 15 cpDNA haplotypes: 
one comprising all sequences for E. madsenii and E. falci- 
sepalum F. Lehm. & Kraenzl., and the other grouping E. 
rhopalostele and most putative hybrids in a single haplo- 
type (Figure 1A, Additional file 1: Table SI). The highest 
haplotype diversity was observed in E. madsenii, which 
had nine haplotypes, whereas the lowest diversity was 
shown by E. rhopalostele, which only had one haplotype 
(H6). The predominant haplotype in E. madsenii (HI) was 
also shared with E. falcisepalum. Three additional haplo- 
types were observed in E. falcisepalum. The majority of 
hybrids shared the haplotype H6 with E. rhopalostele 
(83%), whereas the remaining hybrids contained a unique 
haplotype (Hll; Figure 1A). 

GenBank accession numbers and sequence statistics 
are given in Additional file 2: Table S2 and Additional 
file 3: Table S3, respectively. Hybrids usually showed 
higher nucleotide diversity than their progenitor species 
(Additional file 3: Table S3). 

Genetic composition of hybrid zones 

Bayesian assignment results obtained with STRUCTURE 
confirmed that the allopatric parental populations were 
exclusively composed of pure genotypes because the 
assignment results were very high with a threshold 
q- value > 0.90 (Figure IB). However, in all hybrid zones, 
more than half of individuals morphologically classified in 
the field as parental species showed intermediate ^-values 
(0.10 < q < 0.90) with STRUCTURE (Figure IB). 

Because genome size (GS) values suggested the in- 
volvement of E. falcisepalum and E. madsenii in the ori- 
gin of hybrids in HZ1 (see below), Bayesian assignment 
tests of this hybrid zone were run with these species 
(instead of E. rhopalostele) and hybrids indeed showed 
intermediate ^-values in this population (Figure IB). 
Surprisingly, E. falcisepalum was indicated to be des- 
cended from backcrosses with E. madsenii (Additional 
file 4: Figure SI), which is consistent with the phylogen- 
etic proximity of the two species. 

The individuals revealed by STRUCTURE to be 
admixed were assigned to one of the different hybrid 
classes with NEWHYBRIDS, although the genetic com- 
position varied across hybrid zones. Individuals morpho- 
logically classified as E. madsenii were predominantly 
assigned as backcrossed hybrids in HZ1 and Fx hybrids in 
HZ2 and HZ3. In contrast, individuals morphologically 
classified as E. rhopalostele were assigned as backcrossed 



hybrids in HZ2 and as Fx hybrids in HZ3, whereas individ- 
uals morphologically classified as E. falcisepalum was 
assigned as a pure parental species in HZ1 (Figure IB). 

Hybrids were predominantly assigned as Fj hybrids, 
although F 2 individuals were also present in all popula- 
tions (Figure IB). 

Genome size 

GS values were significantly different between the three 
pure parental species and the mean values were 2.69 pg 
in E. madsenii, 3.13 pg in E. rhopalostele and 4.08 pg in 
E. falcisepalum (i-3,250 = 19.199, P < 0.001). Analyzing the 
individuals according to the genetic groups assigned by 
NEWHYBRIDS allowed us to determine that the GS of 
¥ l hybrids conformed to expectations (GS parentl + 
GSparent2/2) as they usually showed values intermedi- 
ate to those of the parental species in all hybrid zones 
(Figure 2). Backcrossed individuals were usually grouped 
with the respective backcrossed parental species, and F 2 
hybrids had GS values higher than those of the remaining 
genetic groups (Figure 2). 

Some individuals of E. falcisepalum had a GS of 2.91 pg, 
whereas other individuals showed a value of 5.84 pg, which 
suggested the occurrence of polyploidy in this species 
(Figure 2). For the parental species, no significant dif- 
ferences were observed between the GS values of allopatric 
individuals and those assigned by NEWHYBRIDS as pure 
in the three hybrid zones (t = 0.117, df = 98, P = -0.026 for 
E. madsenii; t = -0.026, df =64, P = 0.979 for E. rhopalos- 
tele; t = 0.296, df = 34, P = 0.769 for E falcisepalum). 

Postmating prezygotic isolation: formation of fruit 

Experimental pollinations between the pure parental in- 
dividuals assigned by NEWHYBRIDS yielded different 
results among the three hybrid zones. In HZ1, interspe- 
cific fruit set was lower than intraspecific fruit set, espe- 
cially when E. madsenii was the ovule donor (Figure 3, 
Table 1). In HZ2 and HZ3, no differences were observed 
between intra- and interspecific crosses when E. rhopa- 
lostele was the ovule donor, whereas interspecific fruit 
set was slightly lower than intraspecific fruit set when E. 
madsenii was the ovule donor (Figure 3, Table 1). No dif- 
ferences were observed between populations in the inter- 
action of population x treatment (Table 1). The strength 
of postmating prezygotic isolation (Rlpostmating/prezygo- 
tic) was therefore very low, but variable among species 
(Table 2). 

Postmating postzygotic isolation: development of 
viable seeds 

Among all species, the percentage of viable seeds was 
78.91% ± 25.43 (mean ± SD). No significant difference was 
observed between intra- and interspecific crosses or be- 
tween populations (Table 1). Consequently, the strength 
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(See figure on previous page.) 

Figure 1 Genetic variation in £ madsenii, E. rhopalostele, E. falcisepalum and their hybrids across three hybrid zones. (A) Statistical 
parsimony network of plastid haplotypes based on sequences from six chloroplast DNA regions (trnL-trnF, rps!6, rpoCI, psbK-psbl, rnatk and rbcl). 
A circle's size is proportional to the haplotype frequency. Small empty circles represent single mutational steps. (B) Posterior probabilities (a) for 
the three hybrid zones of Epidendrum madsenii 'MAD' and £ falcisepalum 'FAL' (HZ1) and E madsenii 'MAD' £ rhopalostele 'RHO' (HZ2 and HZ3) 
analyzed with STRUCTURE and NEWHYBRIDS. Individuals identified in the field, based on morphological characters, are delimited by dashed lines. 
Each vertical bar represents an individual. The proportion of color in each bar represents an individual's assignment probability, according to 
different categories (pure parental species, Ft and F 2 hybrids, and the respective backcrosses). 



of this postulating postzygotic barrier (RIpostzygotic) was 
null in all populations. 

Niche models 

The present bioclimatic niche climate was consistent 
with the currently known distribution of E. madsenii 
and E. rhopalostele (Figure 4 and Additional file 5: 
Figure S2). The area under the curve (AUC) score for 
these models was very high (0.88 for E. madsenii and 0.91 
for E. rhopalostele). The projection of the future bioclima- 
tically suitable area showed a general shift of the current 
predicted range because most southern favorable areas 
will disappear (Figure 4). In addition, a severe reduction in 
the area predicted to be suitable for the occurrence of E. 
madsenii and E. rhopalostele was shown (Figure 4). 

Both species will lose more than half of their current 
distribution area between the present day and 2080 
(from present 56,238 km 2 to future 29,131 km 2 in E. 
madsenii and from 71,787 km 2 to 37,439 km 2 in E. 
rhopalostele). Consequently, the future scenario showed a 
decrease of 43.2% in the minimal predicted area (MPA) 
intersection of E. madsenii and E. rhopalostele (from 
43,260 km 2 to 16,827 km 2 ; Figure 4), which suggested a 
decrease in the geographical area potentially suitable for 
hybridization. Nevertheless, the future scenario reveals the 
presence of suitable conditions for the occurrence of the 
three hybrid zones studied here (Figure 4). 

Discussion 

A complex scenario of reticulate evolution in 
Neotropical orchids 

This study has shown conclusively that a recurrent series 
of interspecific hybridization events and backcrosses 
among closely related taxa has lead to the generation of 
genetic novelty and to a complex pattern of ongoing re- 
ticulate evolution within Epidendrum. Our molecular, 
reproductive and genome size data indicate that three 
different species have given rise to two morphologically 
similar hybrids in different areas of Ecuador. In two of 
the hybrid zones studied (HZ2 and HZ3), hybrids were 
derived from E. madsenii and E. rhopalostele, whereas in 
a third area (HZ1), hybrids were instead derived from E. 
madsenii and E. falcisepalum, of which the latter species 
is also of hybrid origin. The bayesian tests assigned E. 



falcisepalum as a backcrossed hybrid with E. madsenii 
and genome size values are consistent with the presence 
of chromosome duplication as a result of polyploidy. Thus, 
not only different species were involved in the process of 
hybridization but in some cases it was also followed by 
polyploidy. Although our preliminary chromosome counts 
indicate that the three species are diploid (2« = 28 chro- 
mosomes; results not shown), we were unable to count 
the chromosomes of all individuals studied, and the differ- 
ences in genome size observed could represent different 
chromosome numbers. 

Our data indicate that E. rhopalostele was in most 
cases the maternal progenitor of the hybrids, except in 
the population HZ1 where E. madsenii seemed to have 
acted as the ovule donor. The fact that most cpDNA se- 
quences were shared with these species is consistent 
with this hypothesis. The predominance of these two spe- 
cies as the maternal progenitor (in HZ1 for E. madsenii, 
and in HZ2 and HZ3 for E. rhopalostele) is also consistent 
with the results of controlled pollinations because fruit set 
was significantly higher when these species acted as the 
ovule donor. 

Hybridization is therefore highly directional and leads to 
the asymmetric formation of hybrids in nature, as re- 
ported in other studies [41]. Asymmetric hybrid formation 
is not unusual in nature [42], and can be caused by com- 
plex genotype-environment interactions [43] or complete 
cytoplasmic incompatibility [44]. Differences in the flow- 
ering period of co-occurring species might also explain 
the presence of a unidirectional gene flow in hybridizing 
populations [45]. Although preliminary, an ongoing study 
reveals differences in the flowering patterns of these popu- 
lations (Vega, pers. comm.). 

Outcomes of hybridization in Epidendrum 

Although hybridization has been long postulated to be an 
important evolutionary mechanism in Epidendrum, genetic 
confirmation of the presence of hybrids in this genus has 
been reported only in the last 3 years and, although results 
are limited to two examples, hybridization is believed to be 
a frequent phenomenon [46,47]. For example, the use of 
nuclear and plastid microsatellites revealed high potential 
for interspecific gene flow even between species with two 
different ploidy levels, i.e. E. fulgens (2n = 2x = 24) and 
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Figure 2 Genome size obtained with PI for three hybridizing species of Epidendrum: E. madsenii 'MAD', £ rhopalostele 'RHO' and 
E. falcisepalum 'FAL' in three hybrid zones (a-c) and in allopatric populations (d-f) according to the genetic groups detected by 
NEWHYBRIDS (colors as in Figure 1B). Values are expressed in picograms. 
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Figure 3 Mean fruit set in experimental crosses within 
(grey bars) and between (white bars) the studied species of 
Epidendrum. Values represent the mean ± SD (W= 100 plants/cross). 
The first letters indicates the identity of the maternal species: 
FAL = £ falcisepalum; MAD = £ madsenii; RHO = £ rhopalostele. 
Crosses with the same letter do not differ significantly (P > 0.05). 



E. puniceoluteum (2n = 4oc = 52) [39]. Although results sug- 
gested the presence of male sterility in Fi hybrids and 
backcrosses, they also revealed wide ecological amplitude 
of the hybrids [39] , which suggested exogenous selection is 
of low importance in this study case. Furthermore, despite 
the fact that hybrids were triploids (2n = 3x = 38 and 40), 
they were able to backcross and produce viable hybrid line- 
ages [48]. 

The prevalence of hybridization was also reported in a 
study comprising 25 sympatric populations of E. calanthum, 



Table 1 Effects of pollination treatment and population 
on fruit set and seed viability of Epidendrum 



Frui set 




SS 


df 


MS 


F 


P 


Population 


0.1671 


2 


0.083 


0.028 


0.383 


Treatment 


11.021 


3 


3.673 


1.252 


4.604 E-13 


Population x treatment 


0.881 


6 


0.146 


0.050 


0.872 


Error 


23.457 


8 


2.932 






Seed viability 




SS 


df 


MS 


F 


P 


Population 


0.261 


2 


0.130 


0.330 


0.894 


Treatment 


0.335 


3 


0.111 


0.282 


0.970 


Population x treatment 


0.672 


6 


0.112 


0.280 


0.860 


Error 


3.561 


9 


0.395 







Table 2 Strength of postmating prezygotic isolation 
(Rlpostmating/prezygotic) in the different hybrid zones 
(HZ) studied for the three species of Epidendrum: 
E. falcisepalum (FAL), E. madsenii (MAD) and E. 
rhopalostele (RHO) 





FAL 


MAD 


RHO 


HZ1 


0.07 


0.40 




HZ2 




0.30 


0 


HZ3 




0.27 


0 



E. cochlidium and E. schistochilum [49]. These three 
diploid species (2« = 28) frequently hybridize when they 
co-occur, although outcomes depend on the species in- 
volved. F x hybrids were frequent in the sympatric popu- 
lations of E cochlidium and E. schistochilum, whereas in 
the remaining populations backcrosses with E. calanthum 
were predominant when this species co-occurred with E. 
cochlidium, but when it co-occurred with E. schistochilum 
backcrosses with the latter species prevailed [49] . 

In contrast to the previous studies, the presence of 
hybrid genotypes was predominant in all sympatric 
populations and most parental species were no longer 
considered as pure species [49]. The same results were 
observed in the current study, especially in the hybrid 
zones where E. madsenii and E. rhopalostele co-existed. 
A high number of viable seeds were produced in our ex- 
perimental crosses and populations were not exclusively 
composed of V 1 hybrids, which would be the case if the 
hybrids were sterile. The existence of one specific haplo- 
type in some hybrids (Hll; Figure 1) is also consistent 
with the fact that populations are not composed only of 
early-generation hybrid genotypes, which contrasts with 
most studies of hybrid zones between deceptive orchids 
[39,50,51]. Rather, populations were also composed of F 2 
hybrids and most likely of backcrossed hybrids. 

Given that the two hybrids studied here are fertile, 
they can easily form a bridge for gene flow between the 
parental species [52]. The studied species are deceptive 
orchids and pollinators usually do not strongly discrim- 
inate species [39], thus they contribute to introgressive 
hybridization through creation of a hybrid swarm and 
the possible merging of species [53]. Probably many 
other examples of hybrid zones among deceptive orchids 
exist, but molecular studies documenting the genetic 
structure of these hybrid zones are scarce, especially in 
the case of Neotropical orchids (but see [39]). 

Evolution of reproductive barriers in Epidendrum 

The high frequency of hybridization in all populations 
studied suggests the general absence of strong interspe- 
cific reproductive barriers. Although not studied here, 
premating barriers are apparently weak in our sympatric 
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Figure 4 Localities sampled in this study (A); Green dots indicate the 3 hybrid zones studied (HZ) while black dots indicate the 
allopatric populations of E. falcisepalum (FAL), E. rhopalostele (RHO) and £ madsenii (MAD). Predictive ecological model of current 
distribution (B) and the future scenario (C) based on the Maxent algorithm of Epidendrum madsenii (blue) and £ rhopalostele (red). The different 
maps assess similarity of niche models between the two species using the intersection of both MPA species. Predominance of one color 
indicates niche differentiation whereas a dark color indicates overlapping of niche models. 



populations owing to the high frequency of hybrid geno- 
types observed. In accordance with this conclusion, sev- 
eral studies have revealed that Epidendrum is usually 
pollinated by a wide number of species, which do not 
show any pattern of fidelity to a given flower species 
(reviewed in [47]). Hybridization and late-generation hy- 
brids are frequent in all Epidendrum species surveyed, 
even between species with different ploidy levels, and al- 
though some degree of sterility is observed, there is always 
a proportion of fertile hybrids with high fitness in a broad 
amplitude of habitats [39]. However, some caution is 
needed when considering the generally low importance of 
reproductive barriers in Epidendrum because our know- 
ledge of the frequency and outcomes of hybridization is still 
limited, even in the species studied. For instance, we only 
found three hybrid zones even though the species studied 
here could apparently co-exist in more sites (Figure 4). 
Furthermore, no studies have evaluated the role of ex- 
ogenous selection in the establishment of hybrid geno- 
types or its association with the mycorrhizal community. 

The occurrence of polyploidy, which has been previ- 
ously suggested to occur in this genus [38], was confirmed 
in our study. Besides hybridization, polyploidy offers a 
good explanation for the high number of Epidendrum spe- 
cies described. Allopolyploid speciation (duplication of 
chromosomes in hybrids) may occur in these orchids by 
several processes, such as somatic chromosome doubling 
followed by selfing to produce a tetraploid, through fusion 
of two unreduced gametes, or by backcrossing of triploids 
with haploid gametes to form tetraploid progeny [10]. 



Can we deal with taxonomic homogenization and 
blurring of species boundaries? 

The widespread hybridization revealed in this species 
complex raises several taxonomic and conservational 
concerns, because it may compromise the genetic integ- 
rity of the parental species to the point of causing local 
extinctions. This is exacerbated by the finding that fu- 
ture niche models predict a severe decrease in the area 
favorable for the occurrence of the parental species. 

With the exception of the polyploid E. falcisepalum, 
which seems to escape the homogenizing effects of gene 
flow, the remaining parental species can no longer be 
seen as pure-breeding species because they consist of 
predominantly backcrossed individuals (E. madsenii in 
HZ1 and E. rhopalostele in HZ2) or ¥ 1 hybrids (E. mad- 
senii in HZ2 and HZ3, and E. rhopalostele in HZ3). The 
virtual absence of pure-bred individuals in these hybrid 
zones implies that genetic swamping has eliminated the 
parental species and led to the loss of species identity 
[53]. It also highlights the fact that, besides morpho- 
logical characters, other taxonomic tools such as mo- 
lecular markers should be used in order to understand 
species boundaries in orchids. 

Previous studies have demonstrated the existence of 
introgression in plants [54-57], which in some cases has 
led to the homogenization of species characteristics and 
the extinction of parental species in very few genera- 
tions. For instance, Helianthus bolanderi suffered a rapid 
reduction in its population size owing to genetic swamp- 
ing by the common invasive H. annuus [6]. Asymmetric 



Marques et al. BMC Evolutionary Biology 2014, 14:20 
http://www.biomedcentral.com/1471-2148/14/20 



Page 9 of 14 



introgression also occurs between Iris fulva and /. brevi- 
caulis where 72% of the /. brevicaulis -like plants have 
alleles introgressed from I. fulva [58]. 

Understanding if such consequences are derived from 
natural or anthropogenically mediated hybridization is 
important for species management and conservation. 
Allendorf et al. [59] categorized hybridization into six 
types based upon the extent of hybridization and the 
natural or anthropogenic status of the hybridization. 
However, in the present study and most other cases it is 
difficult to know whether hybridization was triggered by 
natural causes or not. In some cases, it can even be ar- 
gued that human interventions simply amplify the po- 
tential for hybridization (e.g., through long-distance 
dispersal) or create environments that are more suitable 
for the new hybrid lineages [60] . 

In Epidendrum, hybridization appears to be mainly the 
result of natural causes and is frequently caused by the 
lack of pollinator specificity and the absence of strong 
reproductive isolation barriers [39]. Therefore, in bio- 
logical groups where hybridization is so frequent and 
creates successful interbreeding lineages, two questions 
arise concerning what to do with the species involved. 
Can we even consider them as "species"? Does it make 
sense to restore the historical genetic structure if threat- 
ened species are involved or should we instead target 
our efforts to the new genetic novelties? 

Based on future projection of their ecological niche, 
and given that conditions suitable for these hybrid zones 
will exist, the hybrid zones can therefore be seen as a 
source of raw material for natural adaptive change 
[61,62]. In this light, the definition of species should bet- 
ter focus on traits that lead to adaptation and conserva- 
tion efforts should be targeted towards evolutionary 
processes that generate taxonomic biodiversity instead 
of preserving the taxonomic entities beyond these pro- 
cesses [63]. Although this strategy cannot be applied to 
all organisms, it certainly provides a good framework for 
determining evolutionarily important units [63,64] that 
are worthy of protection and management in species 
with complex reticulate scenarios, such as those in the 
present investigation. 

Conclusions 

Although hybridization is usually invoked as a driving 
force for the high levels of morphological variation in 
Epidendrum, the existence of hybrid zones and the levels 
of genomic admixture remain unknown in most cases. 
Here, we demonstrate that hybridization can be frequent 
whenever two species co-occur, leading to a complex 
scenario of reticulate evolution and the blurring of spe- 
cies boundaries. Given that a high number of viable and 
fertile hybrid seeds were produced, hybridization gener- 
ated genetic novelties in all hybrid zones. The high 



frequency of hybrids and the fact that the parental spe- 
cies were mostly assigned as backcrossed hybrids sug- 
gests low efficiency of reproductive barriers, at least in 
these populations. This conclusion is consistent with the 
general lack of pre- and postmating barriers reported in 
this genus and, therefore, the large proportion of back- 
crosses observed might also be a feature of other 
Neotropical hybrid zones of Epidendrum. The high fre- 
quency of hidden hybrid genotypes described in this 
study suggests that more genetic studies are needed to 
understand the evolution of the genus Epidendrum. 

Methods 

Plant sampling 

Epidendrum madsenii and E. rhopalostele are two epi- 
phytic orchids patchily distributed in the Andean trop- 
ical rainforest. Epidendrum madsenii is endemic to 
Ecuador and occurs in several disjunct localities on the 
eastern slope of the Andes, and E. rhopalostele is near- 
endemic, as it occurs along the Andean Cordillera from 
Central Ecuador to its border with Peru [36]. Flowering is 
reported to occur from November to May in £ madsenii 
and from May to September in E. rhopalostele [36]. The 
sampling was designed to cover the whole geographic dis- 
tribution of these species of Epidendrum (Additional file 6: 
Table S4; Figure 4). All reported hybrid zones were 
sampled (HZ1, HZ2, and HZ3: Additional file 5: 
Figure S2) and, despite intensive fieldwork, no add- 
itional hybrid zones were found. In addition, E. falcise- 
palum found in one of the populations (HZ1) was also 
sampled because this species could also be considered 
as a putative progenitor of the hybrids either by flow- 
ering phenology or morphological traits. Three allopatric 
populations of each species (MAD1, MAD2, and MAD3 
in the case of E. madsenii, and RHOl, RH02, and RH03 
in the case of E. rhopalostele) were also included in the 
study as reference populations (except E. falcisepalum for 
which only one allopatric population was found: FAL1; 
Additional file 5: Figure S2). In total, 286 individuals were 
tagged and sampled: 100 morphologically assigned as E. 
madsenii (70 from the hybrids zones - HZ1:25, HZ2:22, 
HZ3:23 - and 30 from allopatric populations - MAD1:10, 
MAD2:10, MAD3:10), 95 morphologically assigned as E 
rhopalostele (65 from the hybrids zones - HZ1:21, HZ2:24, 
HZ3:20 - and 30 from allopatric populations - MAD:10, 
MAD2:10, MAD3:10), 55 putative hybrids (HZ1:20, HZ2:19, 
HZ3:16) and 36 morphologically assigned as E. falcisepalum 
(26 from the hybrids zones HZ1 and 10 from the allo- 
patric population FAL1; Additional file 6: Table S4). 
Given that these species can propagate by vegetative 
reproduction, individuals were collected with a minimum 
sampling distance of 10 m. Voucher specimens were 
deposited at the herbarium of the Universidad Tecnica 
Particular de Loja. Field studies were conducted in 
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accordance with local legislation. Samples were preserved 
as silica-gel dried leaves and stored at -80°C until use for 
DNA extraction. 

Molecular methods 

DNA isolation, AFLP procedure and sequencing 

For the 286 individuals sampled, DNA was extracted using 
the DNeasy™ Plant Minikit (Qiagen, Hilden, Germany) fol- 
lowing the manufacturer's instructions and stored at -20°C. 

Procedures for AFLP analysis basically followed the 
protocol developed by [65]. An initial screening of 32 
combinations of selective primers was performed using 
five individuals from different populations. Three com- 
binations were selected because they yielded clear and evenly 
distributed bands: EcoRI-AC (FAM)/MseI-CTA, EcoRI-AG 
(FAM)/MseI-CTT and EcoRI-AGG(VIC)/MseI-CTC. 

A reproducibility test was performed by re-extracting 
DNA from two plants per species and per population 
and repeating the whole AFLP procedure. The error rate 
was calculated for every primer combination as the ratio 
of mismatches (scoring of 0 vs 1) over phenotypic com- 
parisons in AFLP profiles and subsequently averaged 
over the three combinations. Non-reproducible frag- 
ments were excluded from the analyses. 

PCR amplifications were performed for 127 individuals 
using an Eppendorf thermocycler in 50 uL reactions con- 
taining 10 ng DNA, lx PCR buffer (Invitrogen, Sao Paulo, 
Brazil), 2.5 mM MgCl 2 , 0.5 ul each primer (10 uM) and 
1 ul Taq (45-130 ug/ml) using 5-10 individuals per popu- 
lation/species. Six chloroplast regions were analyzed: trnL- 
trnF (primers c and f: [66]), rpsl6 (primers rpsl6F and R: 
[67]); and four primers from the CBOL Plant Working 
Group [68]: rpoCl (primers 2f and 4f), psbK-psbl (primers 
F and R), matk (primers 3 F and 1R) and rbcl (primers F 
and R). Primer sequences and PCR conditions were ob- 
tained from the literature (see references above). All PCR 
products were purified using the UltraClean™ PCR 
Clean-up™ Kit (MoBio Laboratories, Inc., Carlsbad, CA, 
USA) in accordance with the manufacturer's protocol. 
Amplifications were performed using the original amplifi- 
cation primers. Purified PCR products were sequenced 
in both directions on a 3730 DNA analyzer (Applied 
Biosystems, Foster City, CA, USA). The sequencing re- 
actions were performed in a total volume of 10 uL contain- 
ing 30-50 ng DNA, 5 uM each primer, 2 uL ABI PRISM 
BigDye Terminator v3.1 cycle sequencing ready reaction 
kit (Applied Biosystems), and 1 uL of 5x Sequencing Buffer 
(Applied Biosystems). 

Data analysis 

AFLP data Amplified bands were aligned with the internal 
size standard using the ABI PRISM GeneScan Analysis 
Software version 3.1 (Applied Biosystems). Subsequently, 
fragments of each primer combination were scored 



automatically with Genographer version 1.6.0 (Montana 
State University, http://hordeum.oscs.montana.edu/geno- 
grapher/) either as present (1) or absent (0), and manually 
corrected. Peaks were recorded in a range from 50 to 
500 bp. Highly reproducible AFLP patterns were found 
for all replicates and the average error rate was estimated 
as 4.2% across all three primer combinations in agreement 
with previous studies below 5% [69] . Eight unreliable frag- 
ments were removed leading to a total of 223 bands. No 
identical multilocus phenotypes were observed between 
individuals or populations. 

To infer population structure and the genetic compos- 
ition of hybrid zones, two Bayesian clustering methods 
were performed using the software STRUCTURE 2.2. 
[70] and NEWHYBRIDS version 1.1 beta [71]. Each hy- 
brid zone was analyzed separately using allopatric popu- 
lations of each species as reference samples of pure 
individuals. STRUCTURE was first used to classify indi- 
viduals as either pure parental species or hybrids using a 
threshold of q > 0.90 to assign pure individuals and 
0.10 < q < 0.90 to classify hybrids, where q represents the 
admixture proportion. We assigned two groups (K = 2) 
in this study because we assumed that two species con- 
tributed to the genetic pool of the sample. Calculations 
were performed running 10 simulations under the admix- 
ture model. The same genotype data used in the STRUC- 
TURE analysis was analyzed with NEWHYBRIDS to 
classify the genetic composition of the hybrids in six dif- 
ferent classes: Fi, F 2 , backcross to each parental species, 
and pure parental species, using a threshold of 0.75. A 
burn in of 50,000 steps followed by run lengths of 300,000 
were used in each program. Similarity coefficients between 
runs and the average matrix of ancestry membership were 
calculated using CLUMPP version 1.1 [72] and visualized 
using DISTRUCT software [73]. 

Sequence data Sequence alignment was performed 
manually using BioEdit 7.0.0 Sequence Alignment Editor 
[74], which was also used to check electropherograms. 
DnaSP version 3 [75] was used to characterize DNA 
polymorphism. Within-species diversity was estimated 
with Nei's haplotype diversity (Hd) and in terms of 
weighted sequence divergence with nucleotide diversity 
(tt) [75]. Chloroplast DNA sequences (the four regions 
concatenated) were analyzed using statistical parsimony 
as implemented in TCS 1.21 with gaps coded as missing 
data [76]. 

Genome size 

Genome size was measured for the 286 individuals that 
were genotyped. Fresh leaves of Epidendrum were co- 
chopped with a fresh leaf of Pisum sativum (internal ref- 
erence standard with 2C = 8.76 pg) using a sharp razor 
blade in a Petri dish containing 1 ml WPB (0.5 mM 
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spermine.4HCl, 30 mM sodium citrate.3H 2 0, 20 mM 
MOPS, 80 mM KC1, 20 mM NaCl, 0.5% (v/v) Triton X- 
100, pH adjusted to 7.0). The nuclear suspension was re- 
covered and filtered through a 50 um nylon filter to re- 
move cell fragments and large debris. Nuclei were stained 
with 50 mg/ml propidium iodide (Fluka, Buchs, 
Switzerland), and 50 mg/ml RNase (Sigma, St Louis, MO, 
USA) was added to the nuclear suspension to prevent 
staining of double-stranded RNA. Five minutes after stain- 
ing, the relative fluorescence intensity of at least 3,000 nu- 
clei was analyzed in a Partec CyFlow Space flow 
cytometer (Partec GmbH., Miinster, Germany), equipped 
with a green solid-state laser for PI excitation, using the 
FloMax software (Partec GmbH). As a quality control, 
only when CV values of G0/G1 peaks were below 3% were 
the analyses saved, otherwise sample preparation was re- 
peated. The mean and standard deviation of the mean 
(SD) of each sample were calculated. The normality of 
the distribution of genome size (GS) of all samples was 
assessed using the Kolmogorov-Smirnov test. Differ- 
ences in GS between hybrid classes were evaluated 
using analysis of variance (ANOVA). In those cases in 
which ANOVA revealed significant differences, the 
Tukey HSD post-hoc test was performed. 

Assessment of reproductive isolation 

Controlled pollinations were performed during January 
2012 in the three hybrid zones and using the individuals 
considered "pure" species by STRUCTURE and NEW- 
HYBRIDS (e.g., q > 0.90). Pollinations were also per- 
formed in one allopatric population per species (MAD1 
for E. madsenii and RHO for E. rhopalostele). Crosses 
were performed between individuals of the same species 
(intraspecific crossings) and between individuals of dif- 
ferent species (interspecific crossings). Pollinations were 
performed by removing pollinia with a plastic toothpick 
and placing them on the stigmas of other individuals 
from the same or the other species. Interspecific crosses 
were performed in both reciprocal directions. To pre- 
vent contaminatory pollination, plants were bagged with 
1-mm-mesh nylon tulle prior to flowering. A total of 50 
randomly selected flowers were used in each treatment. 
For each inflorescence, two flowers per treatment were 
used to avoid potential negative effects of over-pollination 
on fruit set and seed viability. Flowers were monitored for 
fruit set after anthesis and bagged with 1-mm-mesh nylon 
tulle to avoid the loss or predation of fruits and seeds. Fol- 
lowing [39], samples of 300 seeds per fruit were stained 
with tetrazolium and analyzed to quantify seed viability. 
The percentage of viable seeds was determined by dividing 
the number of viable embryos by the total number of em- 
bryos analyzed. Previous studies showed that the percent- 
age of seed viability estimated by the tetrazolium test was 
very similar to the percentage of seed germination in 



several species of orchids [77,78]. The effects of treat- 
ments on fruit set and seed viability were tested with a 
GLM with pollination treatment and populations as fixed 
factors and individuals as a random-effect factor. Fruit set 
and seed set were respectively log- and square-root trans- 
formed before the analyses and, when necessary, analyses 
were followed by ad-hoc comparisons using the Scheffe 
test or the f-test. 

Two measures of postmating isolation and a combined 
index were calculated as follows [79]. First, we estimated 
postmating prezygotic isolation as the proportion of in- 
terspecific fruits formed after interspecific crosses rela- 
tive to the proportion of fruits formed after intraspecific 
crosses: 

Postp rezygotic = l-(fruit set formed in interspecific crosses)/ 
(average fruit set formed in intraspecific 
parental crosses) . 

Likewise, postmating postzygotic isolation was calcu- 
lated as the proportion of viable seeds obtained from in- 
terspecific crosses relative to the proportion of viable 
seeds obtained from intraspecific crosses: 

Postp OS t Z ygotic = of viable seeds formed in interspecific 
crosses) /(% of viable seeds formed in 
intraspecific parental crosses) . 

The strength of postmating barriers varied between 0 
(no isolation) and 1 (complete isolation). When interspe- 
cific crosses performed better than intraspecific crosses, 
the barrier was set to zero. 

Niche models 

Niche models were generated using the Maxent algo- 
rithm [80], which estimates a target probability distribu- 
tion by finding the probability distribution of maximum 
entropy subject to a set of constraints that represent the 
incomplete information about the target distribution 
[81]. The Maxent algorithm was chosen because it is 
specially recommended when a low number of presences 
(<30) is available [82-85]. A total of 36 localities obtained 
from fieldwork were considered: 19 from E. madsenii and 
17 from E. rhopalostele (Additional file 5: Figure S2). No 
records from distribution databases were included because 
of the difficulty of checking the accuracy of the identifica- 
tions. No niche models were constructed for E. falcisepa- 
lum owing to the low number of known localities. GIS 
layers used in the construction of niche models included 
the 19 climate data variables from Worldclim (http:// 
www.worldclim.org/) with a 30 arc-second resolution 
(approximately 1 km 2 ) and altitude obtained from SRMT 
(http://srtm.csi.cgiar.org/) at the same spatial resolution. 
Only uncorrelated environmental variables were consid- 
ered. To select final variables, a jackknife was carried out 
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in order to select variables higher than 0.75 using the 
AUC predicted performance. 

The species' niche was also projected on a future cli- 
mate scenario for the end of the twenty-first century, ac- 
cording to the SRES-A1B, with a scenario reflecting an 
important increase in C0 2 concentration as a result of 
non-restricted CO2 emission (Intergovernmental Panel 
on Climate Change, Special Report on Emission Scenar- 
ios, scenario Al, [85]). Future areas of climatic suitability 
were compared with present-day suitable areas and 
present distribution by calculating the percentage of area 
lost under a scenario of unlimited dispersal. Fitness of the 
models was assessed using AUC of a receiver-operating 
characteristics plot [86,87]. To obtain a binary map, the 
threshold of the resulting models was determined as the 
minimum predicted value that included all records of each 
species, so that a 0% omission error (proportion of ob- 
served presences incorrectly predicted) was attained [81]. 
The resulting maps illustrated the minimal predicted area 
(MPA) predicted for 100% of the presences [88]. By 
assigning a 0% omission error, we ensured that all known 
populations were included in the predicted area. 
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